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The understanding of the physics of laser- 
matter interactions in the strongly relativistic 
regime is of fundamental importance. In this ar- 
ticle, a new mechanism of fast electron genera- 
tion at the vacuum-solid boundary of intense laser 
pulse interaction with overdense plasma is de- 
scribed. It is one that has no analogue in classical, 
non-relativistic laser-plasma interactions. Here, 
conclusive proof is provided that the key contri- 
bution to the fast electron generation is given 
by the zero points of the vector potential. We 
demonstrate that the new mechanism leads to 
scalings for the fast electron energy, which explic- 
itly depend on the plasma density, thus providing 
a new insight into relativistic laser-matter inter- 
action. Furthermore, it is shown that this new 
mechanism provides the dominant contribution 
to the interaction by the injection of energy into 
the overdense plasma delivered by attosecond- 
duration electron bunches. This new understand- 
ing will allow the future generation of a single 
ultra-bright attosecond X-ray pulse by suitable 
control of the laser pulse polarization [ij-Q- This 
process will also allow single pulse attosecond 
electron bunches to be generated that can be 
further accelerated in laser wakefield accelerators 
[8]- [10]. Other applications that will benefit from 
this new insight include laser driven ion acceler- 
ators [11]-[15], fast ignition inertial confinement 
fusion [16]- [18] as well as fundamental studies at 
the intensity frontier [19]- [20]. 

Possible mechanisms for generation of fast electrons 
have been extensively investigated in recent years. A 
well known and widely quoted scaling for the energy of 
the fast electrons was proposed by Wilks [22| ■ This scal- 
ing was observed in several experiments (e.g. 23[ and 
24j), however, other experimental and numerical works 
present a more complex picture (e.g. pjj). In terms of 
the absorption fraction, a recent review of many experi- 
mental and theoretical studies of relativistic laser absorp- 
tion |26j showed that absorption fractions from a few per- 
cent to nearly 100 percent have been measured or calcu- 
lated by numerical simulations. This suggests that there 
are several regimes of absorption. Despite the absence 
of agreement in the laser-plasma community on the best 
scaling or theoretical description of absorption, there is 



a consensus that the ponderomotive jxB-force plays an 
important role in generating relativistic electrons. 

The mechanism described in this article is, by nature, 
a non-ponderomotive mechanism of absorption during in- 
teraction of ultra-intense laser pulses with sharp density 
gradients. This theory predicts the generation of fast 
electrons propagating into the plasma bulk and coher- 
ent x-rays in the reflected from the plasma radiation in 
a single self-consistent framework. The most peculiar 
and novel feature of this theory is that the generation of 
bunches of high energy electrons depends on the existence 
of zeroes in the vector potential and not on fluctuations 
in the incident laser intensity alone. For this reason we 
refer to this new mechanism as the Zero Vector Poten- 
tial mechanism. All previously proposed mechanisms - 
including jxB-heating 27] , vacuum heating 28| , Brunei 
heating [29}, anharmonic resonance [3C| etc - make no 
such claim. This Article describes the new absorption 
mechanism, demonstrates the new scalings for the fast 
electron energy and presents numerical simulations which 
demonstrate the predictions of the zero vector potential 
theory for efficient control of the energy absorption. 

The existence of at least two regimes of absorption dur- 
ing intense laser-plasma interaction can be understood on 
the basis of simple physical reasoning. If we denote the 
frequency of the incident laser radiation as uj then the 
jxB-force oscillates with frequency 2u, which is much 
lower than the frequency of relativistic electron oscilla- 
tions in the overdense plasma (w p = \f 4Trn e e 2 /aom e , 
where n e is the plasma density) . Consequently, the ques- 
tion arises, in which density region does the main in- 
teraction between the incident radiation and the target 
take place? It is clear that if the region where the laser 
frequency is of the order of the frequency of relativistic 
electron oscillations is broad enough, then the interac- 
tion is localised in this region and the jxB-force indeed 
drives the electrons to high energies, which they then 
transport inside the plasma, since plasma effects are not 
fast enough to counteract this ponderomotive force. 

However, a completely different regime takes place if 
the density gradient in the vicinity of the critical density 
is steep. In that case, the pressure of the incident ra- 
diation can significantly compress the plasma and shift 
the interaction area to plasma densities which are much 
higher than the relativistic critical density (Fig.Q]a). As 
a result, in the case of a steep density gradient, the laser- 
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plasma interaction takes place at plasma densities well 
above the relativistic critical plasma density. If this is the 
case, the frequency of the jxB-force is much lower than 
the frequency of relativistic electron oscillations at the 
area determining the energy transfer from the incident 
laser pulse to the plasma. Therefore the plasma is able 
to respond fast enough to compensate the pondcromo- 
tive force by the Coulomb attraction between ions and 
electrons. In other words, the plasma adjusts its state 
adiabatically to this force. 
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FIG. 1: a) For plasmas of steep density gradient the pressure 
of the incident radiation can significantly compress the plasma 
thus shifting the interaction area to hig h plasma densities as 
particle-in-cell simulations 15], 31], 32] clearly demonstrate. 

b) During the adiabatic phase the laser light pressure and 
the electrostatic Coulomb force are in equilibrium building a 
potential well in which the electron oscillations are confined; 

c) As the vector potential goes through zero the adiabatic 
equilibrium is disrupted allowing fast electrons to leave the 
potential well. 

Consequently we are in an adiabatic regime in which 
electrons are not directly accelerated by the jxB-force. 
Instead, the skin layer is in equilibrium between the pon- 
deromotive pressure and the electrostatic Coulomb force. 
These two forces create a potential well in which the elec- 
trons are confined to perform finite oscillations. The po- 
tential well can be visualised if we consider the Hamilto- 
nian of a single electron in the plasma skin layer under 



the action of the field of the laser pulse. Due to conserva- 
tion of momentum in direction tangential to the plasma 
we have 
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where p T is the tangential electron and A T is the tangen- 
tial vector potential. The Hamiltonian of a single electron 
can now be written as 
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revealing that the potential well in which the electrons 
are confined when the laser light pressure and the elec- 
trostatic Coulomb force are in equilibrium consists of the 
electrostatic part <!> vanishing inside the quasi-neutral 
plasma and the linear part connected with the vector 
potential A T outside the plasma slab (Fig. [T]b). 

The equilibrium is disrupted for short times when the 
vector potential goes through zero. Since at these mo- 
ments the ponderomotive pressure vanishes, one of the 
walls of the potential well is suppressed (Fig. Q]c), caus- 
ing a violent reconstruction of the plasma skin layer. 

This reconstruction of the plasma skin layer has 
two important consequences. Firstly, the electrostatic 
Coulomb energy, stored in the skin layer over the adi- 
abatic phase, is freed during the reconstruction in the 
form of fast electrons leaving the plasma in the specu- 
lar direction and the violent electron dynamics leads to 
bursts of coherent x-ray radiation in vacuum in the same 
direction. Secondly, upon returning the fast electrons en- 
ter the plasma, thus transporting energy irreversibly into 
the target. 

It is interesting to note that the process of reconstruc- 
tion can be understood in rather simple terms, owing 
to a fascinating phenomenon which occurs only in rela- 
tivistic skin-layers. Indeed, it turns out that in the re- 
gion which contributes to the interaction of the overdense 
plasma with the relativistic radiation the zero point of the 
vector potential moves with velocity c [33| . In order to 
understand this phenomenon let us consider closely the 
relativistic overdense plasma slab and a laser electromag- 
netic field normally incident onto the plasma for the sake 
of simplicity. Since the slab is overdense it expels the 
electromagnetic field and since the electrons are ultra- 
relativistic the electromagnetic field is expelled with the 
speed of light. On the other hand the transverse mo- 
tion of the electrons is connected to the vector potential 
according to Eq. ([]}. The longitudinal velocity of the 
electrons is related to the vector potential by 
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This means that at the zero of the vector potential 
the electron momentum p vanishes and the relativistic 
electrons which are in the vicinity of this point move 
with velocity c either in the direction out of or into the 
plasma. Thus the ultra-relativistic electrons stick to the 
zero of the vector potential, leading to the zero being 
"dressed" by co- moving relativistic electrons during its 
motion through the plasma. In its turn this overdense 
slab of ultra-relativistic electrons expels the electromag- 
netic field with velocity c, thus closing the self-consistent 
picture of the motion of the vector potential zero. No 
confusion has arisen here between formulating the pon- 
deromotive force in either the E, B or cf>, A formulation. 
One arrives at the same conclusions with either formula- 
tion. 

Note that the longer the distance that the zero travels 
in the skin layer, the higher the number of electrons that 
stick to this zero. The momentum of the trapped elec- 
tron bunch increases due to the electrostatic attraction 
to the ions until the electrons leave the ion background. 
After that the Coulomb attraction to the immobile ions 
starts decelerating the electrons to return them back into 
the plasma. These returning electrons have already con- 
sumed the electrostatic energy of the skin layer. As a 
result they form a short electron pulse of relativistic elec- 
trons injecting energy into the plasma slab. 

It must be emphasized that it is this motion of the 
zero in the vector potential that defines this mechanism 
as having no non-relativistic classical analogue. The rea- 
son for this is that only in a relativistic plasma can one 
have a vector potential with a moving zero in the skin 
layer. In a non-relativistic plasma, the vector potential 
can only exponentially decay in the skin layer. Thus the 
zero vector potential mechanism is a unique product of 
relativistic electrodynamics. 

Note that due to their relativistic motion the coher- 
ently moving electrons attached to the vector potential 
zero radiate in phase, leading to the generation of a large 
number of coherent x-rays during the laser-plasma inter- 
action. These coherent x-rays were observed in a num- 
ber of numerical simulations [l| and experiments 0, H| • It 
was shown analytically by the theory of relativistic spikes 
that the spectrum of these x-rays is universal and obeys 
a power-law decay, where the intensity of the frequency 
ui scales as 0, Q 



In the adiabatic phase of the interaction the electro- 
static response of the plasma balances the radiation pres- 
sure. Noting this, let us estimate the distance Ax by 
which the electron fluid is displaced into the plasma (see 
Fig. [T] a). When the electron fluid is pushed into the 
plasma a non-compensated ion charge Q — en e Ax is 
left at the plasma boundary. The attractive force act- 
ing on the electron fluid from the side of the ions is 
therefore Fi oc e 2 n 2 (Ax) 2 . In the adiabatic phase, this 
force is compensated by the force acting on the electrons 

E 2 

from the side of the laser light pressure F p — '°" r oc 



(aQCLip/e 2 ) (m e c) , thus leading to the following estima- 
tion for the electron fluid displacement 
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Here A is the laser wavelength and S is a dimensionless 
combination of the plasma density n e , the critical density 
n cr and the dimensionless vector potential ao, such that 
S = n e /n cr ao [Hj]. 

The above estimation essentially defines how much en- 
ergy is stored in the 'plasma capacitor' during the adi- 
abatic phase. During the A T = phase the electrons 
discharge the capacitor and get the energy stored in it. 
Thus the energy of the fast electrons can be estimated as 
the energy which they gain when crossing the capacitor 
of width Ax: 



W = eEi as 



-Ax oc e 2 n e (Ax) 2 oc -^(m e c z ). 



(6) 



Eq. ([5]) gives important insight into the absorption of en- 
ergy in relativistic laser-matter interaction. It shows that 
when the interaction takes place at the critical surface 
(n e = aon cr , corresponding to S ~ 1), the fast electron 
energy is proportional to a . This corresponds to the 
ponderomotive regime discussed above, in which Eq. ([6]) 
reproduces the Wilks' scaling. In the adiabatic regime 
however, the interaction takes place at much higher den- 
sities (5^1) leading to 
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up to a roll-off frequency on the order of uJ ro ii oc I Q 
in agreement with experimental observations. 

From these arguments, one can derive a scaling relation 
for the energy of the fast electrons. In what follows we 
will pursue a simple derivation for the case of a step-like 
interface. This can be derived far more rigorously using 
relativistic similarity theory 33j, |3J] ■ 



This explicit dependence of the fast electron energy 
on the plasma density underlines the important differ- 
ence between the zero vector potential mechanism and 
the ponderomotive jxB-mechanism. The density depen- 
dence predicted by the new zero vector potential mecha- 
nism can be demonstrated by a PIC simulation. Fig. ([2]) 
shows that the inverse-linear dependence of the fast elec- 
tron energy on the density can be observed for a broad 
range of plasma densities, thus confirming the predictions 
of the new absorption mechanism. 
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FIG. 2: Log-log plot of the inverse-linear dependence of the 
fast electron energy from the plasma density (Eq. ©), ob- 
served in PIC simulations for a wide range of densities. 



by P oc a 2 + a 2 oc (a 2 — a 2 ) cos 2 (wi) + a 2 the ampli- 
tude of laser pressure oscillations is the same for both 
polarizations. Note that although the light pressure os- 
cillations are identical, the linearly polarized pulse has 
several vector potential zeroes, while the elliptical pulse 
has no zero crossings. 
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We can also estimate the number of fast electrons gen- 
erated, as the number of non-compensated ion charges in 
the plasma capacitor N = n e Ax oc ao\n cr . Note that 
the number of fast electrons generated by this process 
does not depend on the plasma density, but only on the 
intensity and frequency of the laser. 

The physical arguments presented above suggest that 
the zero points of the vector potential are responsible 
for the generation of fast electrons in relativistic laser- 
overdense plasma interaction and thus for the absorption 
of energy in these plasmas. In what follows we apply 
numerical simulations in order to demonstrate how the 
production of fast electrons and coherent x-rays is gov- 
erned by the zero points of the vector potential. For this 
purpose we consider the interaction of a steep plasma 
slab with relativistically intense laser pulses of different 
polarization thus varying the number of zero points of the 
vector potential while keeping the oscillation behaviour 
of the ponderomotive pressure unchanged. 

We first consider the interaction of two laser pulses 
of linear and elliptical polarization with the plasma slab 
of density n e — 80n cr under normal incidence using a 
1D3P particle-in-cell code [TBI, 31 1, 32 1. The wavelength 
of both pulses is X—l^im and their duration is r=5fs. The 
vector potential for the linear polarization propagates in 
vacuum as a plain wave 



a y = eiyoexp (—(x — ct) 2 /2c 2 T 2 ) cos (27r(a; — ct)/X) , 

where a y o — 20, while the components of the laser pulse 
with elliptical polarization are 



a y = a y o exp (— (x — ct) 2 /2c 2 T 2 ) cos (27r(:z; — ct)/X) 
o-z = o-zo CX P ( — (x — ci) 2 /2c 2 r 2 ) sin (2ir(x — ct)/X) , 

where a y o = y/425 and a z o = 5. Note that this choice of 
amplitudes insures that for the laser light pressure given 



FIG. 3: Logarithmic colour plot snapshots of electron plasma 
density as a function of longitudinal electron momentum p x 
and coordinate x for the case of a) linear polarization and 
b) elliptical polarization. Although the amplitude of laser 
pressure oscillations is the same in both cases the fast electron 
shocks disappear in the elliptical case due to the lack of zeroes 
in the vector potential. Colour scale is in units of kgm 2 s _1 . 

Fig.[3]shows the electron distribution function in colour 
as a function of the electron momentum in longitudinal 
direction p x and the longitudinal coordinate x. One can 
clearly see that during the interaction of the plasma slab 
with the laser pulse of linear polarization electron shocks 
are generated in the skin layer and propagate into the 
bulk of the plasma slab twice per laser period (distance 
A/2 apart). At the head of the plasma one can clearly 
observe electrons with positive and negative momenta of 
equal size, representing the oscillations of the electrons in 
the potential well produced by the laser vector potential 
and the electrostatic Coulomb force. 
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FIG. 4: Spectra of coherent x-ray radiation in the signal re- 
flected from the plasma for linear (black line) and elliptical 
polarization (red line). Due to the absence of fast electrons 
in the case of elliptical polarization (compare Fig. [3] b) the 
production of coherent x-rays is significantly reduced. 



Although the oscillations of the ponderomotive poten- 
tial have not changed, when we move to the interaction 
of the plasma slab with the elliptically polarized pulse 
the picture changes completely. Since there are no zero 
points of the vector potential, no electron shocks are gen- 
erated in the plasma skin layer and hence no shocks prop- 
agate into the bulk of the plasma. The elimination of the 
electron shocks in the simulation with the elliptically po- 
larized pulse cannot be explained by the jxB-mechanism. 
Note that one can clearly see a flux of low energy elec- 
trons leaking from the imperfections of the potential well. 
Although both snapshots are taken at the same moment 
of time the leaking electrons have lower velocity and thus 
cannot penetrate into the plasma bulk as deep as the elec- 
tron shocks. 

The importance of these two simulations can be sum- 
marized as follows : We are comparing two cases. One 
where the laser pulse has a strongly oscillating pondero- 
motive force and zeroes in the vector potential, and one 
where we have the same amplitude of oscillation in the 
ponderomotive force, but no zeroes in the vector poten- 
tial. We see that the highly energetic electron shocks are 
only produced when we have zeroes in the vector poten- 
tial. This numerically demonstrates the essence of the 
zero vector potential mechanism, and this cannot be ex- 
plained by the 'jxB' mechanism. 

Note that according to the physical picture presented 
above, the electron shocks are the consequence of elec- 
trons which were able to escape the potential well during 
the times the vector potential goes through zero. Be- 
fore entering the plasma and propagating into the bulk 
in the form of electron shocks these electrons first leave 
the plasma skin layer and radiate coherent x-rays in spec- 
ular direction. Fig. @] shows the spectrum of the reflected 
electromagnetic field. For the case of linear polarisation 




FIG. 5: a) Vector potential of mixed polarization allowing 
only one zero point; b)-d) Logarithmic colour plot of single 
fast electron bunch generated by a single zero of the vector 
potential a), b) Generation of a sharp electron bunch inside 
the plasma skin layer; c) and d) propagation of the bunch 
into the plasma bulk. The relativistic electron bunch expe- 
riences instability during the propagation into the cold bulk 
plasma, e) Zoom of the single electron bunch in b) shows that 
electrons of 2MeV are confined within a time interval of 5as, 
and electrons between 0.5MeV and 2.5Mev are propagating 
within 130as. Colour scale is in units of kgm 2 s _1 . 



this spectrum follows the harmonic decay law (|4]) derived 
in 0] ■ However the x-ray intensity drops dramatically in 
the case of elliptical polarization, as a result of the miss- 
ing fast electrons. 

The experience of controlling the generation of elec- 
tron shocks by means of changing the number of zero 
points in the vector potential suggests an exciting appli- 
cation of this mechanism of control of the laser plasma 
interaction. Indeed, the zero vector potential mechanism 
suggests that if the vector potential of the pulse vanishes 
only once, an isolated electron bunch can be generated. 

In order to demonstrate this numerically we consider 
the interaction of the plasma slab with a laser pulse of 
variable polarization, constructed in such a way as to 
assure that the vector potential goes through zero only 
once. The vector potential components of this pulse de- 
pend on time as 
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hj = a yO exp (— (x — ct) 2 /2c 2 t 2 ) cos (2tt(x — ct)/X) 

i z = a^oexp (-(x - ct) 2 /2c 2 T 2 ) cos (2n(x ~ ct)/X c + n/8) . 



where A c = 800nm. The vector potential is presented in 
Fig. Ha). 

As a result of the interaction of this pulse with the 
plasma slab a single electron bunch is generated in the 
plasma skin layer (Fig. [5] b) and propagates into the 
plasma bulk (Fig. [5]c). Note that the electron bunch is 
moving with ultra-relativistic velocity through the bulk 
of cold plasma, leading to instability of the bunch struc- 
ture (Fig.[5]d). Fig.[5]e) zooms the isolated electron shock 
presented in Fig. Ob) and reveals the spatial width and 
the geometry of this shock. One can observe that elec- 
trons of different energies are localized in an astonishingly 
narrow area, thus all electrons of 2MeV energy are inside 
a 15Aarea (5as duration). This width is much shorter 
than both the plasma length and the Debye length corre- 
sponding to the temperature used in the simulation. This 
is a direct indication to an unusual mechanism of electron 
bunch generation which cannot be reduced to standard 
plasma physics phenomena, taking place on longer scale- 
lengths. In the shock observed the lower energy electrons 
run ahead of the higher energy electrons and that electron 
separation results in a locally monochromatic spectrum 
of the fast electrons, as explained by the zero vector po- 
tential mechanism. Note that when the vector potential 
zero accompanied by fast relativistic electrons leaves the 
plasma there is no spatial separation over energies, and 
relativistic electrons of all energies propagate together at 
about the speed of light. However, the Coulomb attrac- 
tion to the ions turns back the velocities of the electrons 
with lower energies earlier than of those with higher en- 
ergies. As a result, in the electron shock returning to 
the plasma, the lower energy electrons run ahead. In the 
simulation presented on Fig. [5] e) the distance between 
0.5MeV and 2.5Mev electrons is about 130as. The gen- 
eration of a single electron bunch by such an engineered 
pulse cannot be explained by the jxB-mechanism, nor 
can the jxB-mechanism explain the incredibly short du- 
ration of the electron bunches. 

Finally let us emphasize that the zero vector poten- 
tial mechanism applies not only for the case of normal 
but also for oblique laser incidence onto the plasma slab. 
In this case the vector potential zero propagates out of 
the plasma in specular direction leading to a burst of 
fast electrons and coherent x-rays in the same direction. 
Upon their return to the plasma the fast electrons prop- 
agate inwards in the direction of laser pulse propagation. 
This behaviour can be clearly seen in a particle-in-cell 
simulation of oblique laser-plasma interaction (Fig. [6]). If 
one considers the same interaction from the jxB perspec- 
tive then one would argue that the electron bunches must 
enter normal to the plasma surface (the direction of the 
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FIG. 6: Logarithmic colour plot of p x — p y phase space from 
boosted ld3v PIC simulation. Colour scale is in units of 
kgm 2 s _1 . 



gradient in the ponderomotive potential). Only the zero 
vector potential mechanism explains why they propagate 
along the direction of laser pulse propagation. 

In conclusion, we presented a novel mechanism of en- 
ergy absorption in relativistic overdense plasmas and de- 
rived a new scaling for the energy of fast electrons which 
explicitly shows dependence on the plasma density. The 
new non-ponderomotive mechanism explains the relation 
of relativistic absorption to attosecond time scales. We 
demonstrated that the generation of fast electrons and 
coherent x-rays in these plasmas is controlled by the zero 
points of the vector potential allowing the generation of 
sharp shocks of fast electrons propagating into the plasma 
bulk twice per laser period. Controlling the laser pulse 
polarization in such a way that it goes through zero only 
once leads to the generation of an isolated fast electron 
bunch. 
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